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o : 

£\j i We present an exact field theoretical representation of the statistical mechanics of simple classical 

liquids with short-ranged pairwise additive interactions. The action of the field theory is obtained by 
performing a Hubbard-Stratonovich transformation of the configurational Boltzmann factor. The 
mean field and Gaussian approximations of the theory are derived and applications to the liquid- 
. vapour transition considered. 
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I. INTRODUCTION 



This article on a somehow new presentation of the theory of liquids is dedicated to D. Levesque at the occasion of 
his 65 th birthday. Since Dominique and I share the same office, it was nearly impossible to keep secret the writing 
of the paper and I could not help myself discussing with him some parts of the manuscript. His deep knowledge of 
the theory of liquids made his remarks always pertinent and the final status of this article owes very much to our 
. discussions. 

The purpose of this work is to try to apply some of the techniques of statistical field theory to the physics of 
simple liquids. By "simple" we mean a liquid made of a single, chemically inert component which can be satisfactorily 
Ch ' modelled by pairwise additive and spherically symmetric interactions 1 . More precisely we shall consider the case 
of a pair potential v(r) which can be written as the sum of a hard core interaction plus a short-ranged tail —w(r) 
which is supposed to decay faster than l/r 3+e at large r so that the thermodynamic limit exists 1 . Long-ranged 
electrostatic interactions are thus excluded of the present study as well as molecular liquids or mixtures. With the 
help of a Kac-Siegert-Stratonovich-Hubbard-Edwards (KSSHE) 2-6 transformation it is formally possible to reexpress 
the grand partition function of this system of decorated hard spheres as the grand partition function of bare hard 
spheres in an external stochastic potential or "field" <f> with a Gaussian measure 3, . We can thus describe the fluid 
with the help of a field theory characterised by an action Ti.((f>) which however is quite complicated since, besides a 
simple quadratic term, it also includes the grand potential of the hard spheres which is a non local functional of the 
' field (p. Quite immodestly we however assume that this functional is perfectly known. Of course, it would be more 
fN| I satisfactory to treat the hard sphere interactions at the same level as the tail w(r) but the KSSHE transform requires 
a well behaved interaction (namely its Fourier transform must exist and be well behaved, see e.g. the appendix) 
and, unfortunately, hard core or more realistic less singular repulsive short-ranged potentials do not fulfill these 
mathematical requirements. 

All the sophisticated techniques developped in statistical field theory 11-14 can thus be applied a priori to the 
case of liquids. Feynman graphs are no more difficult to compute than Mayer graphs which makes possible to 
perform effectively perturbative expansions of various thermodynamic quantities attached to the liquid in the KSSHE 
representation. Recently this type of approach was used to rediscover the low fugacity 9,10 and high temperature 10 
expansions of the free energy and the equation of state of the restricted primitive model of electrolytes (RPM), results 
which had been obtained years ago before by awkward graph resummation techniques 15 ' 16 . 

Another powcrfull technique of statistical field theory is the so-called loop expansion which amounts to expand 
functionally the action H,{4>) around a saddle-point 12-14 . The lowest order (zero loop) of the expansion, defines the 
mean field (MF) level of the theory. The MF theory of liquids in the KSSHE representation has been studied recently 
by Brilliantov et al. both for homogeneous liquids 7 and for liquids inhomogeneous in one direction of space 8 . In the 
former case, it yields to a van de Waals like theory of the liquid-vapour transition and, in the latter one, it allows 
the calculation of the MF surface tension. In the present paper we reconsider and extend Brilliantov' s work on the 
MF level of the KSSHE theory and go a step further by considering also the one-loop approximation of the theory or, 
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more precisely, a simplified version of it, called the Gaussian approximation in the literature 11 . In the homogeneous 
case it is shown to be strictly equivalent to the random phase approximation (RPA) of the theory of liquids 1 . The 
present work is hopefully a first step forwards a renormalized theory of the critical points of liquids at the one-loop 
order. 

Our paper is organised as follows. In section II we expose in details some important properties of the KSSHE 
transform and establish notably the relation between the correlation functions of the fields and that of the fluctuations 
of densities. In section III the general MF equations of the KSSHE field theory are derived and their properties studied 
in details. As a subproduct of this formalism we obtain exact bounds for the free energy and the grand potential 
functionals. In next section IV the Gaussian approximation is derived and shown to be equivalent the RPA theory. 
In section V we reexpress the action H((f>) under a form similar to that of the Landau-Ginzburg (LG) Hamiltonian 
of magnetic systems which allows us to map the liquid-gas transition on the ferromagnetic one. In this way one 
can identify clearly the order parameter of the liquid-vapour transition which is, at least at the MF level, nothing 
but the density of the fluid, which rules out the field-mixing hypothesis 17 ' 18 . In section VI the MF and Gaussian 
approximations are used to study in details the critical point of the fluid. We conclude in section VII. 

II. THE KSSHE TRANSFORM 

A. The KSSHE Transform of The Boltzmann Factor 

We consider the case of a simple three dimensional fluid made of N identical hard spheres of diameter a with 
additional pair interactions denoted for convenience —w(i,j). Since w(i,j) is an arbitrary function in the core , i.e. 
for rij < a, it is always possible to introduce the following decomposition : 

w(i,j) = w_(i,j) - w+(i,j) , 

w±(k) = J d 3 r cxp(ifc.r) io±(r) >0 (Vfc) . (2.1) 

In order to define safely the KSSHE transformation it is crucial to assume that the Fourier transforms of the pair 
potentials w±(r) do exist and are both positive functions of k (for more details see the appendix). Moreover we 
demand w±(k = 0) and w±(r = 0) to be both well-behaved (finite) quantities. Clearly w+(r) is the repulsive part of 
the tail and W-(r) its attractive part. 

Obviously the decomposition (2.1) is not unique, with the unpleasant consequence that approximate theories can 
depend upon the pecular choice made for w±(r) inside the core. We must precise however what we mean by an 
" approximate theory" . If we mean a systematic expansion of some physical quantity in powers of some small physical 
parameter then the "theory" should be independant of the decomposition (2.1). This point was scrupulously considered 
in ref. 10 where it was shown that each term of the high-temperature and low activity expansions of the equation of 
state of the RPM obtained in the framework of the KSSHE formalism were indeed independant of the regularisation 
of the Coulomb potential inside the core. 

By contrast, the small parameter of the loop-expansion considered in the present work is quite mysterious (we could 
have noted it h according to the tradition !) and certainly not related to any physical parameter of the system 12-14 . 
As a consequence, each term of the expansion of the free energy of the liquid given in the following sections depends 
explicitly upon the regularisation of w(r) at short distances. We shall explain later on how to deal with, and even 
take advantage of the explicit dependance of the theory upon the decomposition (2.1) inside the core. 

We denote by f2 the domain of volume V occupied by the molecules of the fluid and by T their temperature. 
It is convenient to assume that O is a cube with periodic boundary conditions (PBC). In a given configuration 
cj = (N; f\ . . . rjv) the microscopic density of particles reads 

N 

p(r» = £<5( 3 )(r-rl), (2.2) 
»=i 

and the Boltzmann factor can thus be written 

exp (-f3V H) = exp (-/3V HS M) exp(-Nu s ) exp Q (p\MP^j , (2-3) 

where (3 = 1/fcgT (fcs Boltzmann factor) and v$ = /3w(0)/2 is minus the self-energy of the particles. From our 
hypothesis on w(r), vs is a finite quantity which depends on the regularisation of the potential in the core. In eq. 
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(2.3) exp(— (3Vhs(w)) denotes the hard sphere contribution to the Boltzmann factor and we have introduced the 
symbolic notation 

(p\w\p) = / d 3 f 1 d 3 f 2 p(fi\u)w{ri,r 2 )p{r 2 \u) . (2.4) 
Jn 

With our assumptions the two quadratic forms (p\w±\p) are both positive definite and we can take advantage of the 
properties of Gaussian functional integrals to rewrite 



cxp (p\/3w-\p)^ = (cxp((p|^_))) (3M ,_ 



where the scalar products (p\f±) of the microscopic density pwith the two real scalar fields (p±(f) are defined as 

@\<p±) = [ d 3 f p{r\w)<p±{r) . (2.6) 
Jn 

The brackets (. . -)p w denote Gaussian averages over the fields <p±(r) i.e. 

Af 0w± = J IV±(r>xp(-i(^ ± |(/? W± rV±)) , (2.7) 

where the inverse of (3w± must be understood in the sense of operators, i.e. 

/ d 3 r 3 w±{r u r 3 )w± 1 {r 3 ,r 2 )=5 {3) {r u r 2 ) . (2.8) 
Jq 

One could object that the relations (2.5), (2.7) and (2.8) are rather formal but in simple cases - i.e. if PBC are assumed 
- explicit and unambiguous expressions of the measure T>ip±(r) can be given. See appendix A for more details. In 
order to simplify the notations further, we introduce the complex field ip(r) = tp-{r) + iip+if), define the measure 

Xy(f) ee ZV+(r)XV_(r) , (2.9) 

and the Gaussian average 

(..V=^/^-"^(-^E(^i(^) _1 i^)) ' ( 2 - 10 ) 

where Np w = Np w+ x Np w _ . We can thus rewrite the Boltzman factor (2.3) as 

exp(-/3VM) = cxp(-/3Vff S H)exp(A^ 5 ) (exp((^|p))) /3w . (2.11) 

Admittedly these notations can be misleading but they can be fully justified, see the beginning of section V for a 
more detailed discussion. Eq. (2.11) defines the KSSHE transform of the Botzmann factor. 

B. The Physical meaning of the auxilliary fields 

We have seen in section II A that, in a given configuration u, the Boltzmann factor of the fluid can be written as 
the statistical field partition function 

z(w) = exp(i < p\(3w\p>)=M^ J ZVexp(-%]) , (2.12) 

where h[ip] is a Gaussian Hamiltonian which reads 
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%] = ^E(Ve|(/9«; e )- 1 |^)-( V _|p)-i( V+ |p) . (2.13) 



2 ^ ± 



The saddle point of h[tp] is denned by the set of equations 



Sh 



5ip-(r) 

the solution of which is 



Sh 



= (2.14) 



<P+(ri) =i0 d 3 r 2 w+(fi,f 2 ) p(f 2 ) , 
Jn 



= iPw + (n,.)*p(.) , 
jp_{n) =p W -(n,.)*p{.) , (2.15) 

where the symbol V denotes a convolution in space. Note that eqs (2.15) imply the compact formula 

Tp(n)=pw(f u .)*p(.) , (2.16) 
where Tp(r) — Tp_(f) + iTp + (f), or, more explicitely 

N 

Tp(r) = J2[3w(f-fi) . (2.17) 

i=l 

Therefore, at the saddle point, the KSSHE field Tp(r) is simply minus the local potential energy at point f. Moreover 
one notes that 

hffi = -\ (p\Pw\p) . (2.18) 

i.e. the value of h at the saddle point h\0\ coincides with the configurational energy. Let us make the change of 
variables tp± (r) = Tp ± (r) + Sip± (r) , one has obviously, as a consequence of the stationarity, 

%] = h[Jp] + 1 -Y J (SVeKPwJ- 1 ^) (2.19) 

The positivity of both operators w^ 1 confirms that (Tp, h\0\) is indeed a saddle point. For a more complicated 
Hamiltonian than h[tp], the approximation consisting in truncating the corrections to the saddle point value h[tp] at 
the Gaussian level is called the Gaussian approximation 11,12 . This approximation is obviouly exact for z(u) because 
h[ip] is a quadratic form. Indeed a direct calculation shows that 

z G ee AT-' J VS^(r) exp (-hffl (^K^)" 1 !^)^ = cxp (-hffl) ee z(lu) . (2.20) 

Note that it follows from eq. (2.15) that Tp + is imaginary whereas we assumed the reality of the field ip + . The integrals 
on 8(p + in cq. (2.20) correspond therefore to a mere translation for the contour of integration. 

C. The KSSHE Transform of The Grand Partition Function 

Henceforward we shall work in the grand canonical (GC) ensemble. We denote by p, the chemical potential an 
by ip(f) an external potential. The local chemical potential will be noted v(r) = [3{p — ip(r)) and the GC partition 
function reads : 

oo i r N 

^\ = Y.m dl --- dN e M~PV(uj)) H cxp (v (i)) , (2.21) 

AT=0 ' ^ n i=l 

where i = n and di ee d 3 fi. S can be symbolically rewritten as 
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E[i/] = J dn{co)exp(-(3V(co)) exp((z/|p)) . (2.22) 

With these notations averages will be obtained by the formula 

, A( \\ _ f^M^Mexp(-^H) exp(Hp)) . . 

As well known, S[i/] is log-convex functional of v 19 ~ 21 . Inserting the expression (2.11) in eq. (2.21) one obtains readily 
a result which seems to have been established for the first time by Siegert 3 , i.e. 

S[i/] = (E HS [V]) 0w . (2.24) 

where V(r) = v(r) — vg + f(r). Eq (2.24) tells us that the grand partition function S of a fluid of decorated hard 
spheres is equal to the mean value of the grand partition function Ehs of bare hard spheres in the presence of an 
external stochastic field with a Gaussian weight. 

To make some contact with statistical field theory we also introduce the effective Hamiltonian (or action) 

H[<P] = \ E (<Pe\(PWe)- 1 \<Pe) ~ ^HS M , (2.25) 

which allows us to write alternatively for 5 

SM = J V<p(r) exp(-W[ V ]) . (2.26) 

It will be important in the sequel to distinguish carefully two types of statistical field averages, the already defined 
< . . . >p w and the < . . . >u that we define as 

< >W = JVW) exp(-WM) ' (2 - 2?) 
With these definitions in mind one notes that for an arbitrary functional of field A[~p\ one has 

< A <p >«= — — -r- . (2.28) 

( ^HS W\) f3w 



D. Preliminary Results on Correlation Functions. 

The ordinary and connected correlation functions of the fluid will be defined in this paper as 1,22,23 

c<">M(i,...,n) = = L £ E M. . , 

=l[v\ 01/(1) . . . ou(n) 

= (f[m\^) , 

\l=l / GC 

<- 9 » 

Our notation emphasizes the fact that the (connected and not connected) are functionals of the local chemical 
potential v(r) and functions of the coordinates (1, . . . , n) = (fi, . . . , f n ). We know from the theory of liquids that 22,23 

Gi n) M(l n)=G(")H(l n)-J]n G< m > [v)(h, . . . , i m ) , (2.30) 

where the sum of products is carried out over all possible partitions of the set (1, . . . ,n) into subsets of cardinal 
m < n. In standard textbooks 1 the n-body correlations are usually defined as functional derivatives with respect to 
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the activity rather than with respect to the chemical potential. It yields to differences involving delta functions. For 
instance for n = 2 and for a homogeneous system one has 

G^{v}(l,2)=p 2 g(r 12 )+p5(l,2), 

G^[i^](l,2) = p 2 h(r 12 ) + pS(l,2) , (2.31) 

where p is the density of the fluid and g(r) the usual pair distribution function at chemical potential v, finally h = g — 1. 

Of course we would like to relate the correlation functions of the microscopic density p with the correlation functions 
of the KSSHE field ip defined as 

G^M(1 n) = (n^)) , 

\i=l / H 

G%l[v]{l,...,n) = Gl[v]{l,...,n)-Y J R G<£>[i/](ii, . . . ,i m ) , (2.32) 

From the expression (2.26) of S and the definition (2.29) of G^ one infers that 

G (n) M(l> • • • ,n) = S- 1 M£ J V<p cxp (-1 <^I(«-V>) 2 ffs (77) G%1[V}(1, ...,„), (2.33) 

which yields, with the help of eq. (2.28), to the simple relation 

G(")H(l,...,n) = (G^[F](l,...,n)\ . (2.34) 



We are only half the way since the hard sphere correlation functions G^g [F] are of course complicated functionals of 
if (through their dependence upon V). Let us first examine the simple case n = 1 ( G^ l=1 ^ = p). We first rewrite eq. 
(2.34) for n = 1 as 

p[v]{l) = { P Hs[v]{l)) n , 

= S- 1 ^ 1 /^exp^-l2<^l^«)-V«>) (2-35) 



Then one remarks that 



^_(1) ^+(1) ' ^- dbJ 

which allows us to replace the functional derivative of with respect to i/ in the r.h.s of eq. (2.35) by a derivative 
either with respect to tp- or with respect to ip + . Let us do it in details in the former case. Since for functional 
integrals one has 12 

^||=0, (2.37) 
An integration by parts yields the results 

*](D = -H- as; /*v SmW , 

= (0t»_)-'(l,.). (»>_(•))« . (2.38) 

Repeating the above derivation with the field <p> + instead of <p>- and making use of the compact notations previously 
defined one finally arrives at the set of relations 

<¥>-(!)>« = ^-(1. ■) * (PhsM(.)) h > (2-39a) 
<¥>+(!)>« = i/Ju>+(l, •) * (p fl ^](.))„ , (2.39b) 
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(<p(l)) n = Pw(l,.)*(pHsM{.)) 



(2.39c) 



For a homogeneous system one therefore has (<p) w = 0w(O)p. It can be noticed that eq. (2.39c) is analogous to a 
relation derived by Callen for an array of spins Si 2 , namely < Si >=< tanh(/?/ij + flJijSj) > (hi magnetic field on the 
site i, Jik coupling constant between sites i and j). Callen's exact relation can be used to obtain the MF equations of 
spin systems by neglecting all correlations, which yields the well-known MF equations Si = tanh(/3/ij + (3JijSj) where 
Si is the MF magnetisation of site i. Proceeding in an analogous way with eq (2.39c) on obtains the MF equations 
for the KSSHE field ip: 

W)=Pw(l,.)* P HsW-v s +m.) , (2.40) 

which will be derived more rigorously in next section. 

This method of functional integration by parts can be used a priori at all orders so that to derive relations between 
the correlations of the field p and the correlations G^ of the microscopic density p. We just quote the result 
for n = 2. 

Gf H(i, 2) = /Mi, 2) + /Mi, •) * g {2) MO. •) * /M. 2 ) . 

2) = /Ml, 2) + /Ml, •) * ^ H(, .) * /M, 2) . (2.41) 

We leave the derivation of cqs. (2.41) by the method of integration by parts as an exercise for the reader. The result 
will be rederived and generalised to the case n > 3 by a more elegant method in section V. 



III. MEAN FIELD THEORY 

A. The grand-canonical free energy 

We define the MF theory or saddle point approximation by the equation 

E MF (v)=exp(-H(<p)), (3.1) 

where at tp = Tp the action H is minimum. It there are several local minima then one retains the absolute minimum. 
The stationary condition reads 



sn 



= (e = ±) . (3.2) 



5ip e (r) 

It follows readily from the expression (2.25) of TL that these eqs. may be written 

Tp_(l) = (3w-(l,.) *p H sp]{-) , 

Jp + (l)=i0w + (l,.)* P HsM(.) ■ (3-3) 

By combining linearly cqs. (3.3) one recovers the MF eq. (2.40) for ip. The similarity between the MF eqs. (3.3) and 
(2.40) in one hand, and that which give the saddle point of h[cj>] in a given configuration of the phase space (cf. eqs. 
(2.15)) on the other hand, is striking although expected. 
Making use of eqs (3.3) it is then easy to show that 

logS M FH = logS H5 [77] - 1 (pHSp]\ftw\pHsW\) ■ (3.4) 

The MF density at point fi is obtained by taking the functional derivative of logSM^M with respect to y(l). It 
yields 

5 log SmfM 

PmfM(1) = 



Mi) 

PHsim) + (phsM | A - fhv * . (3.5) 
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The scalar product in the r.h.s of eq. (3.5) vanishes as a consequence of the stationarity conditions (3.2) and, finally, 
one finds 

Pmf[v]{1) = Phs[v]{1) (3.6) 

We are now in position to compute the MF grand-canonical (MFGC) free energy as the Legendre Transform of 
logSM^M with respect to the local chemical potential v(r) 

0Amf[pmf] = (pmfW) - log S M fM (3.7) 
One finds after simple algebra that A[p], as a functional of p, can be expressed as 

0Amf[p\ = l3A HS [p] + I d 3 rv s p(r) - \ (p\0w\p) , (3-8) 
Jn z 

where -4h-s[/°] is the exact GC free energy functional of the hard spheres at the density p. For a homogeneous system 
Amf(p) is merely the MF free energy of the system in the GC ensemble. 



B. The Mean Field Density Functionals 



We recall first some important properties of the functionals logS[^] and A[p\. Under quite general conditions it can 
be shown that the logarithm of the exact grand-partition function log is a convex functional of the local chemical 
potential v(r) and that the exact GC free energy A[p] is a convex functional of the density p(r) 19-21 . It must be 
stressed that, for a finite system, A[p] differs from the canonical free energy -4c [p] which is not a convex functional 
of p(r) notably in the two-phases region due to intcrfacial effects. Moreover logS[f] and A[p] constitute a pair of 
Legendre transforms, an important property which can be expressed as 

pA[p}=sn P (( P W)-logE[u]) (Vpeft), (3.9) 

log3[i/] = sup«p|i/)-MM) (Vi/GW). (3.10) 
pen 

It can also be shown that the sets of physical densities 1Z and that of local chemical potentials U are both convex 
sets 19,21 . The so-called Young inequalities follows directly from eqs (3.9), (3.10) and read 

pA[p] + IogH[i/] > (p\v) (Vp G H,Vv e U). (3.11) 

Finally we recall that a necessary and sufficient condition for the convexity of logS[^] and A[p] is that their second 
order functional derivatives are positive operators, i.e. 19-21 

(H^^(-G?')MI*')>o(*> el O, 

Moreover the two-body direct correlation function [p] is minus the inverse of the connected two-body correlation 
function [v] i.e. 

GM[v}(l,.)*CW[p}(.,2) = -6(1,2) (3.13) 

where it must be stressed that v[p] is the (unique) local chemical potential corresponding to the density profile p. We 
have now at our disposal all the mathematical tools necessary to prove the two following theorems 

• if w+ = (attractive case) PAmf[p] is an upper bound of /?.4[p] 

• if w_ = (repulsive case) (3Amf[p] is a lower bound of (iA[p\ 
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Lot us consider first the attractive case (w+ = 0). We start with the fundamental relation (2.21) which reads in this 
case 

SH = (exp (logSH 5 [^ - u s + ^})) fjw _ (V*v e U) . (3.14) 
We now apply Young's inequalities (3.11) to logS ffs which yields 

H[i/] > (exp (-0A\p] + (p\v -v s + V-)))p w _ (Vi> G U, Vp e H) 

= exp(-0AHs\p] + {pW-v s ) + ^{p\0w-\p)] (V^eW,Vpeft). (3.15) 



2 

where we have made use of the property (2.5) of Gaussian integrals. Taking the logarithm one arrives at 



pA MF [p] = PAhsIp] + (pWs) - \ ip\pw-\p) > (pW) - logSH; ; (Vi/ e W,Vp e K) . (3.16) 



which implies 



I3Amf[p] > sup ((p\u) - IogS[i/]) = MM («>+ = 0, Vp G ft) . (3.17) 



Good approximate functional @Ahs[p] are available in the literature for the hard sphere fluid. Our (3Amf[p\ could 
be of some use to deal with systems of decorated hard spheres. Note that PAmf[p] is also a functional of the pair 
potential @W-{r) in the core (r < a); therefore an optimized version of the MFGC free energy can be obtained by 
minimizing 0Amf[p] with respect to ty_(r) in the core. 

The repulsive case (io_ = 0) is more tricky. One first remarks that, in this case, one has 

3[i/] = (exp {\og~ HS [v -u s + iV+]))p w+ G «) • (3-18) 

which implies 

E H s[v] = (expQogZ[v + t's + <P+]))f lw+ (Vw G W) . (3.19) 
We are thus led back to the previous case. We therefore apply now Young's inequalities (3.11) to logS which yields 

Z H slis}>e X p(-(3Alp} + (p\v + v s ) + l(p\f3w+\pf) (V^W.VpeK), (3.20) 

where we made use, once again of the fundamental property of Gaussian integrals (cf eq. (2.5)). Taking the logarithm 
we are thus led to the following inequalities 

PA[p]>-\ogZ H sW\ + {p\v + vs)-\{p\pw\p) (Vi/GH.Vpeft), (3.21) 

where we have noted that, in the repulsive case, w + = —w. Since the inequality (3.21) is valid for any v G U for a 
given p G 1Z it is also true for v* = sup v (y eM) which yields the desired result 

PA[p\ > (3Amf[ P ] (w- =0, Vp G K) . (3.22) 

In the repulsive case an optimized version of the MFGC free energy can thus be obtained by maximizing /3Amf[p] 
with respect to w+(r) in the core. 

In the general case, i.e. when w = W- — w + (w- ^ and w_ ^ 0) Amf[p] is neither an upper or a lower bound 
for the exact GC free energy. Moreover, in the homogeneous case, there is no clear relation between Amf[p] and the 
Gibbs-Bogoliubov bounds 1 for the free energy 

C. The Mean Field Correlation Functions 

1. The connected two-points correlation function 
The MF connected two-points correlation function will be defined according to eq. (2.29), i.e. as 
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(2) _ 8 2 log~ MF [v] 

G MF >](1,2)_ Sv[l)5v{ 

_ 5 P mfW}(2) 

~ Su(l) ■ [6 6) 

Since P mf[^]{2) = Phs[^ — + ^](2) (cf. eq. (3.6)) the MF density pmf\v\ is a functional of the local chemical 
potential v directly and also through the KSSHE fied Tp[v\. Therefore one has 

C7^, c H(l,2) = G^ iC [F](l,2) + J rf3G^ >c [F](l,3)g||. (3.24) 

We note that it follows from from the MF eq. (2.40) that 

§=«3,.)*4>](.,2) (3.25) 

from which we infer the relation 

G%M = G ( hLM + G% C [V] * 0w * G$ FiC W] (3.26) 
which can be solved formally to give 

G^fM = (1 - * G^JF])" 1 * G^ jC [F] . (3.27) 
For a homogeneous system the Fourier transform of G^ F c therefore reads as 



G%M(k) ggy^j . (3.28) 



2. Tfte direct two-points correlation function 
The MF two-points direct correlation function will be defined according to eq. (3.12), i.e. as 



C MF [p\(l,2)- 5p{1)5p{2) ■ 



In the terminology of statistical field theory C y MF [p] (1, 2) is a vertex function, it is related to the usual direct correlation 
function c(l,2) by the relation Cj^ F [p](l,2) = cmf[p](1,2) — <5(1, 2)/ i o(l) 22 ' 23 . It readily follows from the expression 
(3.8) of the MFGC free energy that we have 

= aP»[p](l,2)+Wl,2) . (3.30) 
Some remarks are at order. Firstly, it can be checked that the mean-field operators —C^ F is indeed the inverse of 

(2) 

Gmfc' more precisely one checks readily that, formally 

C$ F \p](l,.)*G$ FtC [v](.,2) = -5(1,2) , (3.31) 

where it must be stressed that in the above equation p denotes the mean-field density at the local chemical potential 
v, i.e. p = phs\v\. For a homogeneous system the Ornstein-Zernike (OZ) eq. (3.31) can be written simply in Fourier 

~(2) ~ (2) 

space as G^p c [v](k)C MFc [p](k) = — 1. Secondly it must be stressed that the expression (3.27) of the connected 
correlation function G^ Fc is rather formal. Its validity, and thus the validity of the OZ relation (3.31) as well, is 

guaranteed only is the inverse of the operator 1 — (3w * G$ s c [v] does exist. It can be expected to be true only if 
T > T c where T c is the (mean field) critical temperature. Below T c the first derivatives of the MF grand potential 
with respect to v have discontinuities at the chemical potential corresponding to the coexistence of the vapour and 
liquid phases and the second derivatives (i.e. G^ F c [v\) are not defined. 
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3. Rigorous bounds for Log S[v] 

(2) 

It transpires from the above discussion that for T > T c the MF correlation function G MFc [v] is a positive operator, 
i.e. that 

(Vi/.VfceH) (h\G% c [v)\h)>Q. (3.32) 

It follows then from the OZ eq. (3.31) that for T > T c and for all densities p E TZ, C MFc [p] is also a positive operator. 
Consequently, both log SmfM and (3Amf[p\ are convex functional of respectively the generalised chemical potential 
v and the density p. Moreover they constitute a pair of Legendre transforms by construction and thus satisfy to the 
eqs. (3.9), (3.10). Exact bounds for logS[z/], valid for T > T c , can be deduced from this property 

Let us start with the case of a repulsive pair potential (i.e. W- =0). In this case (3A[p\ > {3Amf[p\ for all p as 
shown by eq. (3.22). Therefore, for T > T c 

logE MF [v] = sup ((p\v) - (3A M f[p\) (W E U) , 
pen 

> sup ((p\v) - (3A[p\) (Vi/ E U) , (3.33) 
pen 

from which we conclude that 

logS MF H > logSH (T > T c , w_ = 0, Vi/ E IX) . (3.34) 
In the case of a purely attractive pair potential (i.e. w+ = 0) it can be shown similarly that 

log Smf \y\ < logSH (T > T c , w + = 0, Vis E U) . (3.35) 



IV. THE GAUSSIAN APPROXIMATION 
A. The grand potential 

Let us define the Gaussian approximation of the KSSHE field theory in the following way 11 . First we write 

V- = <P- + £- , 

y+=^++C+, (4.1) 

where ^_ and ^ + are real scalar fields (£ will be defined as ^ = + and then we expand functionally the action 
H[ip] (cf eq. (2.25) ) up to second order in £± around the MF solution. In this way the exact action H[f] is replaced 
by an approximate action Hg [f] defined symbolically as 



H G [v\=H[Tp]+ l -jj{l)d{2) 



S 2 H 



Mi)M2) 



_ C(l)€(2) • (4.2) 



Note that the terms linear in ^ are absent from eq. (4.2) as a consequence of the stationarity condition (3.2). Taking 
into account the explicit expression (2.25) of H we have more precisely 

where the operators A± are functional of V = v — i>s + V and are defined as 

A ± 1 ^(/3 W± )- 1 ±G^ C [77]. (4.4) 

or equivalcntly 

A± =pw±*(l±G ( $ StC [v}*[3w±)- 1 . (4.5) 
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Clearly the Gaussian approximation makes sense only if both operators A± are positive, which will be assumed 
henceforth. 

The Gaussian grand partition function is obtained by replacing Tt[<p] by its Gaussian approximation TLgW\ m the 
expression (2.26) of E[v], which gives 



E G [u] = exp (-Hffl)M£ J V£ + cxp fcl^lk) + i (z + \G$ s jm-)^ 



(4.6) 



S G \v] can be explicitly computed in the case of a homogeneous system. For instance, by first performing the Gaussian 
integrals on the variable then on , one easily obtains that 

E G [v] = exp {-H[Tp]) A/a+ M x _ , (4.7) 
or, by reverting the order of integrations on £ + and , that 

E G M = exp {-HM)N^ A/"a_ . (4.8) 

The operators X± which enters eqs. (4.7,4.8) being defined by the relations 

X ± l = Ag 1 + G%ljv] * A T * . (4.9) 

Since Af/3 W = Np w+ Np w+ we see that only ratios of the normalization constants TV of various operators enter the 
expressions (4.7, 4.7) of S G [i/]. These ratios are easily evaluated with the help the eq. (A12) of appendix A . For 
log S = VfiP holds for a homogeneous system, one obtains finally for the pressure P G 

d 3 k ( A_(k)X+{k) \ 



f}P G {")=fiP MF {v) + \ |^Alog^ 



Pw+{k)/3w-(k) I 

= PPhs(v)-1(3w(0)p 2 hs (v)+ 1 - J ^tog(l-w(k)G™ tC [v](kj) . (4.10) 



B. The correlation functions and the free energy 

The mean values of the type (2.28) with the full action H replaced by its Gaussian approximation TLg will be 
denoted by < ... > G . Clearly one has < £± >g— from which we infer that 

PgW] = PmfM = PhsP] • (4.11) 

The various pair correlations of the KSSHE fields <p± are obtained by applying Wick's theorem (cf eq. (A14)) with 
the result 

(¥>±(1)¥>±(2)) GiC =X ± (1,2), (4.12a) 
( V+ (l) V _(2)> GiC = iA+(l, .) * G$ S JV}(., .) * (^_(.)<p_(2)> G c . (4.12b) 
After some algebraic manipulations more easily effectuated in Fourier space, one obtains the more transparent relations 

<¥>±(l)p±(2)> GiC = ^±(1, 2) T 0w ± {l, .) * gSp iC M(., •) * dw ± {.,2) , (4.13a) 
{<p+{\)<P-{2)) G<c = ipw + (l,.)*G$ F Jv](.,.)*f3w-(.,2), (4.13b) 
(^(1M2)) G C = W (1, 2) + (3w(l, .) * G$ F Jv}(., .) * (3w(., 2) . (4.13c) 

(2) (2) 

The comparison of these formula with eqs. (2.41) shows that G MF = G G a well known property of the Gaussian 
approximation 11 . Recall however that we have defined G^ F as the second functional derivative of Emf M with respect 

(2) 

to v. The second functional derivative of ^g\ v \ with respect to v differs of G G as defined above by terms involving 
the three-body correlation function of the HS reference fluid. This difference between the two possible definitions of 
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(2) .... . (2) 

G G coincides in fact with the so-called one- loop correction to G G which play an important role in the renormalized 
theory of the critical point (such a theory remains however to be worked out explicitly !). The previous remarks could 
be extented to all correlation functions; for instance at the one-loop order the density does not coincides with the MF 
or Gaussian result. 

Since pc = Pmf the Legendre transform of logScH, i.e. the free energy of a homogeneous system in the Gaussian 
approximation is easily derived. Let fa = Aa/V be the specific free energy; we have 

Pfa(p) = PfMF(p) + /3A/(p) , (4.14a) 

PfMF(p) = Pf HS (p) - ^(0)P 2 + ^ PW (0) , (4.14b) 

PAf( P ) = \ J ^0 log (l - (M(k)G%> StC \p]W) • (4.14c) 
It should be noticed that the expressions derived above in the framework of the Gaussian approximation for the 

^(2) (2) 

correlations C G , and G G c , and for the free enncrgy fc coincide with those obtained in the framework of the RPA 
theory of simple fluids 1 . As well known, the RPA theory becomes exact in the limit of an infinitely long-ranged 
interaction w(r) — j 3 wo(jr), 7 — > 1 . 

Note that, as in the case of the MF theory, the free energy (4.14) depends on the value of the pair potential w in the 
core, which is arbitrary. We can take advantage of this freedom of the theory by requiring the stationarity condition 

<)IJ<: - for I |rl I < cr . (4.15) 



5j3w{r) 

For it is easy to derive from (4.14) that 



<Wg p 2 ,~ 



5[3w(r) 2 

the stationarity condition (4.15) on [3 fa implies the nullity of the radial pair distribution function in the core, an 
important physical requirement. This optimized Gaussian approximation coincides of course with the ORPA theory 
of the theory of liquids 1 . 

V. THE LANDAU-GINZBURG ACTION 

We first note that, quite remarkably, all the results of the previous sections could have been obtained more simply 
by ignoring the complications introduced by the presence of two distinct KSSHE fields (p + and associated with 
the repulsive and attractive part of the pair potential respectively, and by considering rather a unique KSSHE field 
ip. Indeed, taking as a starting point the simplified definitions 

SM = Ay J V<p(f) exp(-WM) , (5.1a) 

= \ (vKPwy 1 ^) - , (5-lb) 

Mfh» = J V<p{r) exp(-((^|(/? W )-V)) (5.1c) 

without any requirement on the positivity of the quadratic form < ipKPw)^ 1 ]^ >, yields the MF and Gaussian 
expressions for the MF and Gaussian expressions for log S, the free energy and the correlation functions. This remark 
justifies a posteriori the notations (2.9) adopted at the beginning of section II A. Henceforth, in order to avoid any 
awkward complications in further developments we shall adopt definitively the somehow abusive definitions (5.1). 

We show now that the action Ti[ip] can be rewritten as the Landau- Ginzburg (LG) action of a magnetic system 
in the presence of a magnetic field B. Following Brilliantov 7 we make choice of a reference chemical potential vq for 
the HS reference system (supposed to be uniform for simplicity and to be specified further) and perform a functional 
Taylor expansion of logS^s around vq. Let us first define for further convenience 

Av(r) = v{r) - u s - vq , (5.2a) 
<j)(r) = Av(r) + <p(r) . (5.2b) 
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With these definitions we have v = v$ + <j>. Now, assuming the analiticity of logins at we have, for v in some 
neighborhood of uq 

log Sjts [F] - log 3*3 M + Whs\M) + \ (^GhIMU) - VM , (5.3a) 

oo 

V/M = -E^i /dl---dnG^ iC [^](l,... ) n)^(l)...0(n). (5.3b) 

n=3 

Therefore one can rewrite the expression (2.25) of as 

H[i/] = exp(-W) S LG [S] , (5.4a) 

W= ^(Ai/IOSti;)- 1 ^!/) -logHirsN , (5.4b) 

where Ti is a simple Gaussian functional of t/ and Slg[.B] may be seen as the partition function of a LG model in the 
presence of a magnetic field B. Indeed 

E LG [B] = M£ J V<f>(r) exp(-W LG ,oM + (B\4>)) , (5.5a) 

HlgM = i (0| A" 1 \4>) + Vz[$ , (5.5b) 

where the free propagator of the theory A is given by A -1 = {fiw)^ 1 — G$ s [vq] or equivalently 

A = l3w*(l~l3w*G ( H ) s Jv })- 1 , (5.6) 

and the magnetic field B by 

B = PHS [v \ + (J3w)- 1 *Au . (5.7) 

It can be noticed that the interaction term Vx[4>] of the LG action TClg.o is non local and does not exhibit the usual 
symetry Vx[4>] = Vx[—<f>] of Isink like systems. The correlation functions of this field theory are defined as usual as 

G^[B](l,...,n) = /n^)\ 

\ i=1 I Hlg 



1 S n E LG [B] 
E LG 6B(l)...SB(n) 
6 n \ogE LG [B} 
SB(l)...SB(n) ■ 



(5.8a) 

>)r D ,„ \ 5 n \og3 LG [B] 

G\'[B]{l,...,n) = — — — — . (5.8b) 



Since it follows from the definition (5.7) of B that 

we deduce readily from the relation logShv] = logS#s[i/o]+logSLG[-B]— < Ai>| (flw)^ 1 1 Az/ > /2 that p = (/3w) _1 *(< 
>h L g — Ai>) = ((iw)^ 1 * < <p >Ui which coincides with relation (2.39c). Performing then n successive functional 
derivatives of logS[i/] in taking into account eq. (5.9) one proves easily the following relations 

G^[u] = -{Pw)- 1 +G { Z[B] (5.10a) 

G^M = G ( ;.ilB} forn>3, (5.10b) 

where 

G ( £[B}(1, . . . , n) = f dV . . . dn'iPw)- 1 ^, 1') . . . (0w)-\n, n') Gj»[B](l', . . . ,n') . (5.10c) 
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Since the KSSHE fields and (p differ by a constant Av the connected correlation functions G^ C [B] and G^ c [v] are 
equal for n > 2. Note that in the case n = 2 eq. (5.10) coincides fortunately with eq. (2.41). 

When evaluated for a homogeneous field $ = cte the LG Hamiltonian HhaWi = ^lg[0] where /ilgM takes the 
familiar form 11 ' 14 

oo 

h LG[4>] = 2 r 0</> 2 - B <t> + ^2 ' ( 5 - H ) 

n— 3 

where r = K — K , with if = l/(/3u>(0)) proportional to the temperature and K = Pjj S (vo)- In this case we have 
simply B = phs{vq) + KAv and finally u n = — /?P^(z/o) = — p^g ^(^o)- Note that, in general, odd and even powers 
of 4> enter the expression of ft-LG [</>]• 



VI. THE LIQUID- VAPOUR TRANSITION 



A. The mean field KSSHE theory 



1. General discussion 



We apply the results obtained on the previous sections to discuss now the liquid-vapour transition. The subject 
was discussed to some extent by Brillantov in ref. 7 with however some inaccuracies that we correct and some missing 
points that we include in our discussion. 

With the notations of section V the MF equations (3.2) take the familiar form 

SHlg,o 



= B = p HS [M + {P™y 1 *^- (6-1) 

Eq. (6.1) gives the MF magnetisation of the magnetic system associated to the fluid in the presence of the magnetic 
field B. For a homogeneous system we seek a solution cj> = <fi = cte of eq. (6.1). This is equivalent to find the minima 
of the function h(<j>) = h LG ((f>) + KAv 2 /2 - f3P HS (v ). We have 

h(4>) - yAt, 2 + y^ 2 - KA^ - l3P HS (v + 4>) ; (6.2a) 

ti(<f>) =K(f>-KAi/- P Hs(yo + <t>) , (6.2b) 
h"(<j>)=K-p' HS (is + <f>) . (6.2c) 

Moreover the link between the fluid and the associated magnetic system is made through the relation 

B = phsM+KAv . (6.3) 

Let us first discuss the case K > w(0) > 0, i.e. when the pair potential is more attractive than repulsive. It is 
clear from eqs. (6.2) that the convexity of the function (3Phs{ v ) governs the properties of the function h(<p). Since 
0Phs(v) is convex its second derivative PP'hs(^) = p'hs( v ) 1S positive for all v. Moreover for all reasonable equation 
of state of the HS fluid the function PP'^ s {v) exhibits a unique maximum at some v§ — as can be seen in figure (1) 

where the graphs of some derivatives p^ s (u) of the HS density are sketched within the framework of the Carnahan- 
Starling (CS) approximation 1 . In the CS approximation we have i/q = —0.025 and PP'nsi^o) = = 0.090. Similar 
values are found in the Percus-Yevick (PY) approximation 1 (via the compressibilty route which is simpler, see table 
I) since, at moderate densities, the predictions of the two theories are in close agreement. 
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FIG. 1. The HS density pus(v) and its first derivatives p^ s {v) versus the chemical potential v. The graphs were drawn 
with the help of a parametric representation \p^\{v) = a(rj),v — b(ri)] (rj packing fraction) extracted from the CS equation of 
state 1 . 



Therefore, for K > Kq we have h"{4>) > for all values of <j> which implies that h{4>) is convex. Since \im ( j ) ^± 00 h(<f>) = 
+00 then h{4>) has one minimum for some <f> and this minimum is unique, vq being given, let us choose now v such 
that B = then it follows from eqs. (6.3) and (6.2b) that the MF equation h'(<p) — can be rewritten as 

K(j) = Phs{vo + 4>) - Phs{vo) ■ (6.4) 

A solution of cq. (6.4) is obviously = and it is the unique solution, for K > Kq. For K < Kq the function h(4>) is 
no longer convex and other solutions arise. 

Let us assume now that K < Kq . We examine first the case = Noting for convenience (f>* the KSSHE field 
and h* the Landau function in this case we have from eq. (5.11) 



h*((j)*) = yAi/* 2 + -j-<f>* 2 - B*4>* + -^cj)* 4 + 0(<t>* 5 ) ; (6.5a) 

B* = Phs(vo) + KAv* , (6.5b) 

Au* = v — vs — ^0 1 (6.5c) 

< - -PhsH) > • (6.5d) 
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(2) (3) 

Since = —Ph's( u o) = an< ^ moreover u\ = —PHsi^o) > ( see the table and figure (1)) the Landau function 
h*((f>*) describes a 2 nd order phase transition with a critical point at (B* = 0, K = K ). For K < K and B* = 
the solution <j> = of eq. (6.4) becomes unstable and two stable solutions ±<f) ^ corresponding to a liquid and a 
vapour phase emerge. 

Now if ^ Vq a term in <fi 3 is present in the Landau function h(<p) at B = which describes now a first order 
transition without critical point. This apparent paradox is solved in the next section. A sketch of the functions h*(4>*) 
and h{<j>) is given in figure (2) to illustrate our discussion. 



h(<|>) 




FIG. 2. The Landau-Ginzburg Hamiltonian h(<f>) in zero field (B = 0) for various choices of the reference chemical potential 
vo and for different temperatures K in the vicinity of K c . In each bunch of curves K increases from bottom to top. h(4>) 
was computed in the CS approximation. Only the choice vo = yields a second order phase transition with a critical point 
(bottom curve). In this case the dashed symmetric curves where obtained by truncating h(4>) at the order (f> 4 (included). For 
vo / Vq, h(4>) describes a first order phase transition (top curves) without critical point. The units of h(4>) were not specified 
for more clarity. 

Let us now consider briefly the case K < of a repulsive tail. In this case h^ 2 \(f>) is negative for all <j) an d the 
function h(<j>), ^ e M is concave with a unique maximum for some <fi. For <j) assuming imaginary values h(<j>) is 
obviously a minimum of h(<j>) and the point (<p, h(<j>)) is indeed a saddle point. Therefore for K < no transition 
occurs, at least at the MF level. Henceforth we shall assume that K > 0. 
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2. The Hubbard- Scho field transform and the critical point 



When we make choice of the special value v = v$ the Landau function has the following exact expression 



h* (</>*) = ^-Ai/* 2 + i</>* 2 - KAv*<j>* - (3P HS (u* + 0*) , 



(6.6) 



where we have kept the notations of section (VI A 1). A Taylor expansion of eq. (6.6) around (f>* — gives back the 
usual ferromagnetic LG Hamiltonian (6.5a) with no term in (j)* 3 . The following set of transformations 



obviously enables us to write h(<f>) = h* ((/>*) where h((f>) is given precisely by eq. (6.2a). The relations (6.7) constitute 
the exact Hubbard-Schofield transformation which allows to get rid of the term in <fi 3 in h((f>) 5 . Note that under this 
transformation 



as can be infered readily from the definition (6.3) of the magnetic field B. As a consequence, the function h{4>) at zero 
field which describes a system undergoing a first order phase transition may be seen as the LG Hamiltonian associated 
with a system undergoing a second order phase transition in presence of the field B* = —K(i/q—i>o)—phs{ 1/ o)+Phs{ 1 'o)- 

Henceforth we shall adopt for definitness the choice vo = i/q and drop all the subscripts "*" to simplify the notations. 
With this choice, the liquid is associated with a magnetic system with magnetisation <j> and the liquid- vapour transition 
and its critical point can be made in a one to one correspondancc with the para-ferromagnetic (second order) transition 
and its critical point. Since the magnetisation and the density of the fluid are equal (up to a multiplicative constant) 
the order parameter of the liquid-vapour transition is merely the density. The field-mixing hypothesis of Rchr and 
Mermin 17 ' 18 is therefore meaningless in the present theory. 

The MF critical point of the liquid-vapour transition is therefore defined by the conditions 



The first condition (6.9a) can be satisfied only if w(0) = (w-(0) — w+(0)) > 0, i.e. if the pair potential is more 

attractive than repulsive. K c = p#s(t/o) ~ 0.290 (see the table) is a universal quantity in the sense that it does 
not depend on the pair potential, however the critical temperature T c — K c /w(0) is not universal and depends in 
particular upon the regularisation of w(r) in the core. The second condition (6.9b) determines the critical density 
since it follows from (3.6) that p c — phs{ v o) ~ 0.249 (see the table). p c is therefore also a universal quantity. Finally 

the third condition (6.9c) determines the critical chemical potential Av c = — Phs{ u o) / Pijs^o)- ^- v c is universal but 
not v c which once again depends upon w and its regularisation. 

All the conclusions concerning the universality of some of the critical parameters in the MF approximation originate 
in our choice of treating at the same level the attractive and the repulsive part of the pair interaction. If the repulsive 
part w + is treated exactly and included in the reference system then p c , K c , and Av c are no longer universal quantities. 

In rcf. 7 an approximate Hubbard-Schofield transformation was devised by considering the truncated Taylor expan- 
sion of h{<f>) rather than its exact expression (6.2a); it yields p c ^ pnsi^o) and that p c can moreover be adjusted at 
will to fit with experimental or numerical simulation results; these two conclusions are however incorrect. 



(j) = (j)* + Vq - i/ , 

At/ = At/* + i/q - i/o , 



(6.7a) 
(6.7b) 



B = B* + K(i/q - u Q ) + phsM - Phs{vq) , 



(6.8) 



Kc = K , 
4 = 0, 
B c = 0. 



(6.9a) 
(6.9b) 
(6.9c) 



3. Mean Field solution for T > T c 



We have seen in section (VI A 1) that for K > K c = p$ s (i/o), the MF solution <f> is unique and satifies the relation 



h'(4>) = 0^=>K<p = KAv + phs(vo + <t>) ■ 



(6.10) 



The general results of section (III) can be specialised to the uniform case and yield 
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PPmfM = 0Phs(vo + (f>)- 7^Phs("o + <t>) , (6.11a) 
Pmf{v) = phs{vo + 4>) > (6.11b) 

2 

Pf MF (p) = Pf HS (p) - y/J«KO) + f /M«) • (6.11c) 

Note that the critical pressure 

/3Pc = /3Ph S M--% M - ( 6 - 12 ) 
2p^(^o) 

is a universal quantity (one finds for the critical compressibilty factor Z c = f3P c / p c ~ 0.3589(C5') , 0.3598(PFC)). 
Note also that eq. (6.11b) implies that the critical isochore p = p c is defined by the condition B = 0, since in the case 
the (unique) solution of (6.10) is <fi = 0. 

We complete the above results by establishing now the MF expressions for the internal energy the specific heat 
and the compressibility. The excess internal energy by unit of volume u ex is given by 

U ex (P, = IPfMF(p) - PfHs(p)} , 

= | W (0) - ^-w(0) , (6.13) 

from which it follows that in the MF approximation the excess specific heat vanishes above T c : 

du ex 



a 



V,ex 



0(3 



= . (6.14) 



p 

The isothermal compressibiliy \t = flp~ 2 (dp/ dv)p can be rewritten as 

P (l) , -tn d~4> 

Ph S {v + <t>) dv 

For it follows from the MF equation (6.10) that 

K 



dip 



K - PhU^O + <t>) 



(6.15) 



(6.16) 



we have finally 



XT = l^. ^aki±4 (617) 

P 2 hs( v o + 4>) K-p ( £ s (v Q + <t>) 

Alternatively, \t could have been obtained from the Fourier transform at k — of the two-body correlation function, 
by making use of a well-known expression of the theory of liquids 1 



XT = — G M F,cM( k = °) i 



PhS ( , V + 4>) 

(3 G^sM+W) 



P% s ("o + <1>) l-K-^G {2) Sc [v + 4>]{Q) 



(6.18) 



where we made use of eq. (3.28). Noting that G^ s c [v + 4>](0) = Phs( u o + 4>) we are led back to eq. (6.17). Along 
the critical isochore (p = p c <f> = 0) one has simply 

*M> " | T^k • < 6 - w > 

with a divergence for K — > K c + with a MF critical exponent 7 = 1 as expected. 
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The behavior of Gmf, c (&) f° r k — > along the critical isochore is also easy to study. As a consequence of eqs. (2.31) 
and (3.28) one has 

G M F.c[pc}(k) = z p F P ° - , (6.20) 

where <Ss[p c ](fc) is the Fourier transform of the usual direct correlation function of the HS fluid at the density 
Pc = Phs{vo)- At low fc's one has 

CHs[Pc](k) +0w(k)) =chs[Pc](0) +0w(O)) - ak 2 + 0(k 4 ) , 

a= zl [ ^(2) (0)+ ^ [pc](0)]; (6 .2i) 

which implies, putting all things together and after having noted that 1 — p c CHs\fic](0) = Pel 'Kq that 

G MF M{k) ~ TT 1^ for k - , (6.22) 

where the correlation length £ reads as 

(6.23) 



(if - ifo)V2 



Of course £ is defined only if a > which put some restrictions on the pair potential (cf eq. (6.21) and note 
that c^g[p c ](0) ~ 1.331 (PYC) ). £ diverges for K — > if c + along the critical isochore with, as expected, a MF 
exponent v = 1/2. At the critical temperature eq. (6.22) implies that GMF.c(k) ~ fc~ 2 the usual behavior in the MF 
approximation yielding r\ = for the Fisher exponent. 

We end this section by a discution of the behaviour of the order parameter along the critical isotherm K = K c . 
Near the critical point the MF eq. (VI A 3) at T c can be reexpressed as 

6 = Au+ PHs{Uo+ ^ 



PHsivo) 

1 (3) ( \ 

= A, A„ c + + I ^f^/ + 0{f) , (6.24) 

yielding 

f = —{y-v c )p% ) s {v ). (6.25) 

Therefore along the critical isotherm we have finally 

P = Phs(vo + 0) ~ PhsM + P^(^) 4/3 (6/w 4 ) 1/3 (^ - ^ c ) 1/3 (6.26) 
yielding the classical value S — 3 of the critical exponent, as expected. 

4- Mean Field solution for T < T c 

Below (and near) T c it is sufficient to consider the truncated expansion (6.2) of the Landau function h{<j>). The 
solution of the MF equation is of course well known in this case. Let us define B = B/\B\. Besides <j> = 0, the MF 
equations h^\(f)) — have the solution 11 ' 12 ' 14 



= (j> Q B + 5<f> 
3!\ 1/2 

U4 



S ~^2(K^ + 0iB2h (6 ' 27) 
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In zero field (i.e. B = 0) the solution <j> = is unstable and the two solutions <f> = ±|0 O | are stable. The densities 
corresponding to the magnetisations ±|0 O | are those of the coexisting liquid and vapour; they are given by 



Pi, g = PHsivn ± | </>o |) 



31 x 1/2 



~ PhsM ± P^M (^-J (K K) 1 ' 2 , (6.28) 

yielding a value /3 = 1/2 for the exponent of the order parameter, as expected. It can be noticed that eq. (6.28) 
supports the law of rectiligned diameters with the simple result (valid in the vicinity of T c ) pi + p g = 2 x p c . The 
chemical potential at the coexistence is given by Av coex = — pnsi^o)/ K as a consequence of the condition B — and 
the pressure at coexistence reads as 

pPcoe* = PPhsM - || + ^(K - K ) 2 (6.29) 

Therefore v coex and (3P coex are both regular functions of the temperature K near K c . 

The compressibility is not defined in the two phases region but one can compute the excess specific heat at constant 
volume thanks to a formula due to Yang and Yang 25 . Along the critical isochore one has 

d 2 P d 2 u 

Crri u * coex rr ' J H'coex i a on \ 

V.ex = i —7y^ Pel , (6.30) 



which yields after some algebra to 



ZU4 



The comparison of eqs. (6.14) and (6.31) show that the specific heat at p c has a discontinuity of universal value 
Cv,ex as expected. 



ACy.ex ~ Zp { H ] s (v ) 2 /2u A - 1.0496 (CS) or 1.0074 (PYC), yielding a MF value a = for the critical exponent of 



5. The Gaussian approximation 

As well known, the critical exponents of the Gaussian model are the same as those of the MF theory except a. We 
shall thus content ourselves to compute Cy, e x along the critical isochore above T c . Since Cy. ex {p c ) vanishes in the MF 
approximation above T c , it is equal to (d 2 (3Af /df3 2 ) Pc in the Gaussian approximation where j3Af has been defined 
in eq. (4.14c). Working out the derivatives one find that 



C v , ex (Pc) = ^ I 7^3 G^ F Jp c }(k) (3w(k) . (6.32) 



1 r d 3 k 

2 J Wf 

— (2) 

The integral diverges in the limit K — ► K c + because of the singular behavior of G MF c [p c ] (k) at small k's. The 
method for extracting the divergence of Cv,ex(pc) is well documented 11 ' 12 ' 14 and will not be repeated here. The idea 
is to make the change of variable k' = where £ is the correlation length and to inject in eq. (6.32) the small "k" 

— (2) 

behaviour of G HS c [p c ] (k) as given by eq. (6.22). One finds for K — > K c + 

C vM ~ - *r 1/2 , (6-33) 

yielding a divergence in t~ x l 2 as expected. Below T c a similar behavior (with a different prefactor) can be obtained 
by making use of the Yang- Yang formula. As well known, a comparison of the discontinuity AC'v.ex(pc) of the specific 
heat in the MF approximation (cf eq. (6.31)) and its Gaussian behavior (6.33) allows an estimate of the Ginzburg 
temperature range AK, in our case one has AK ~ a 3 / 2 It depends on the pair potential w through the parameter 
a defined at eq. (6.21). Since a ~ 1 and 114 is tiny in general, AK/K C > 1. 
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VII. CONCLUSION 



In this paper we have presented an exact field theoretical representation of the statistical mechanics of a simple 
model of liquid. The action of the theory is obtained with the help of a KSSHE transform of the Boltzman factor 
of the fluid. Some complications of the formalism such that the fact that the action of the KSSHE theory is not a 
local functional of the field have their origin in the necessity to take an explicit account of the hard core interactions. 
However, it is possible to establish the relations between the correlation functions of the field and that of the density 
of the fluid. Moreover the KSSHE action can be interpretated as the Landau- Ginzburg action of a ferromagnetic 
system in the presence of a magnetic field B which is related linearly with the local chemical potential of the liquid. 

The MF and Gaussian approximation of the KSSHE theory can be worked out explicitly For homogeneous fluids 
the Gaussian approximation coincides with the RPA theory. These two approximations yield a rough description of 
the liquid-vapour transition. The exact Hubbard-Schoficld transformation can be establihcd which yields the usual LG 
Hamiltonian in (p 4 of the theory of critical phenomena. The density emerges as the order parameter of the transition 
which rules out the field mixing hypothesis. The critical density turns out to be independant from the details of the 
attractive part of the potential which is a serious flaw. Moreover several critical properties either do not depend on 
the potential w(r) or depend strongly on the regularisation of w(r) in the core which is in both cases unsatisfactory. 
By lack of place we did not discuss the full one-loop order theory where fresh difficulties arise such as unavoidable 
corrections to the magnetic field B, the order parameter, and more generally all the correlation functions. A fully 
renormalizcd theory at the one-loop order seems however feasible despite these difficulties. 

The MF-KSSHE theory also works for inhomogencous fluids and yields density functionals which could be usefull in 
applications. Moreover for two classes of pair potentials, the MF-KSSHE free energy functional constitutes a rigorous 
bound for the exact free energy functional which could serve as a test in numerical studies. A similar result holds also 
for the grand-potential functional above the critical point. 
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APPENDIX A: FUNCTIONAL INTEGRATION 

In the case where the domain occupied by the particles of the fluid is a cube of side L with periodic boundary 
conditions ( ft = C 3 ) we give an explicit expression of the measure Vip which enters the functional integrals considered 
in this paper. 

In cubico-pcriodical geometries the microscopic density p can be written as a Fourier series 

= JlJ2Pk exp(ifc.r) , 

fceA 

N 

h = Yl cx pH^) > ( ai ) 
i=i 

where A = (2w/L) Z 3 denotes the reciprocal lattice. Note that since p is real one has 

h = PU • (A2) 

The pair potentials w± (rij) between two particles i and j must take into account the interactions between the 
periodical images of i and j. They are both periodical functions of nj which can be written as 1 

w ±(r) = W ±(r + Ln) , (A3) 

n£Z 3 

The above expression for w± (r) can also be written as a Fourier series, the Fourier coefficients of which are precisely 
the Fourier transforms of w± (r) as a consequence of Poisson's formula. Therefore one also has 
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w± 3 (f) = Ta iw±(fc) exp(ifc.r) . (A4) 



L 3 

fceA 



where u>±(fc) is real and has been defined at eq. (2.1). It follows from eqs (Al) and (A3) that, in C 3 , the potential 
energies (j) ^ \w± \p ci ^ in C 3 can be reexpressed as 

1 he 3 \ w f tfO \ 1 ^ w ± (k)\p % \ 2 , (A5) 



2 \ r ' ± lr / L 3 

keA' 

where the sum in the r.h.s runs over only the half A* of all the vectors of the reciprocal lattice A (for instance those 
with n x > 0) as a consequence of the reality of w±(k) and of the symetry relations (A2). 

At this point we introduce a periodic real field <p(r) supposed to be expressible as a Fourier series 



ip(r) = ^3^] tp(k) exp(ifc.f) , (A6) 

and define the measure 14 



L 3 

keA 



Vip=H d 2 !p(k) , 
keA' 

d 2 p(k) = d 5R^(fc) d 3£(jfe) , (A7) 

where it is understood that the domain of integration of both the real and imaginary parts of ip(k) is the whole real 
axis ] — oo, +oo[. 
Defining formally 



keA 

= isE^(W 2 + ^) 2 ) , (A8) 



keA 

and noting that 

®f) = h £ {MP&Wft + ^(Wg) > (A9) 
feeA* 

one deduces easily from the definition (A7) and from properties of ordinary Gaussian integrals 

r + oo , i 

/ dx cxp(--Ar 2 +Ux) = (2ir/A) 1 / 2 exp(--A- 1 J 2 ) (A > 0) , (A10) 
J — oo 2 2 



the fundamental relations 



exp -(pl^l?) = V / . 7 ^ = (exp ((fly))) c , , 

\2\ /; /^cxp(-i(^i(/3^)-i|^)) 

exp (--(pl^l?) = V ) - , 7 ^ = (exp (i (fly))) c, , (All) 

V 2\ /; /^exp(-i(y|(/3^)-i|y)) ^ 

already mentioned in the main text (cf eqs. (2.5)) with less awkward notations. An important result concerns the 
ratio of two Gaussian partition functions. Let Mj be the constant 

Af f = Jvy exp (^ir 1 ^)) , (A12) 
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where / is some positive two-body operator. Then, for / and g positive we have 

/(*) 



N f /M g = n 



~ exp ij/ w io6 i^- < ai3 » 

In practice only functional integrals of Gaussian functionals and their derivatives can be performed explicitly. All the 
feasible integrals can be deduced from Wick's theorem 1 1 ~ 13 which states that 

(<p(l)...<p(2ri)) y =^3/(ii,i2)---2/(i2n-i,i2n) , (A14) 
where the summation runs over all the (2n — 1)!! distinct pairs (ii,^)- 
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TABLE I. Derivatives of the HS pressure at the chemical potential vq where PP^s reaches its maximum. The numerical 
estimates were obtained in the framework of the CS and PYC aproximations. 





CS 


PYC 




-0.025 


-0.059 


PPhs{vq) 


0.433 


0.425 


p HS { Vo ) = (3P^o) 


0.249 


0.246 




0.090 


0.0896 


(2) ( \ 
PhsM 


0.0 


0.0 


4 (3) , n 


0.0116 


0.0119 




-2.766 


-2.743 
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